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Abstract 

7r + 7r~ transitions of heavy quarkonia, especially the Y(3S) — > Y(lS)ir + Tr~ 
decay process, are revisited. In the framework of the Chiral Unitary Theory 
(ChUT), the S wave tttt final state interaction (FSI) is included. It is found 
that when an additional intermediate state with J p = 1 + and / = 1 is intro- 
duced, not only the tttt invariant mass spectrum and the cos9* distribution in 
the T(3S) — > Y(15')7r + 7r~ process can simultaneously be well-explained, but also 
a consistent description for other bottomonia ir + Tr~ transitions can be obtained. 
As a consequence, the mass and the width of the intermediate state are predicted. 
From the quark content analysis, this state should be a bbqq state. 
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1 Introduction 



In recent years, more and more decay data of heavy quarkonia have been accu- 
mulated, and new information on hadron physics has been extracted. Many investi- 
gations along this line, for instance, the decay properties of the ty(2S) — > J/^5/n + n~, 
T(2S) -> T(lS)7r+7r-, T{3S) -> T(2S)n + n~, and T{3S) -> T(1S)ti + 7i- processes 
have been carried out. A commonly used method for such studies is the QCD mul- 
tipole expansion method proposed by Gottfried [2] and further developed by others 
[3, 4, 5, 6, 7, 8]. It was shown that although most data of the mentioned processes 
can be well-reproduced, the nn invariant mass spectrum and the angular distribution 
of the T(3S) — > T(lS)n + n~ decay cannot satisfactorily be described. 

Phenomenological models [9, 10, 11, 12, 13] are also used in such studies. For 
instance, in Ref. [11], Mannel et al. constructed an effective Lagrangian on the basis of 
the chiral perturbation theory (ChPT) and the heavy quark non-relativistic expansion. 
Under the approximation in the limits of the chiral symmetry and the heavy quark mass 
[17], the measured nn invariant mass spectra of the ty(2S) — > J/^!n + n~ , T(3S) — > 
T(2S)n + n~ and T(2S') — > Y(lS)n + n~ decay processes were very well fitted, but not 
of the decay T(3S) — > T(lS')7r + 7r _ . In order to explain the data of dT(T(3S) — > 
T(15 , )7r + 7r~)/(im 7r7r , the nn S-wave FSI was included by using a parameterization. 
However, although the coupling constant ratios (<?2/<?i)&6 i n the T(2S) — > T(lS')7r + 7r _ 
and T(3S) — > Y(2S)n + n~ decays are approximately equal to each other, they are 
much smaller than that in the T(3S) — > Y(lS)n + n~ decay. The later one is about 
10 times larger than the former. It seems unnatural. Moreover, M.-L. Yan et al. [12] 
pointed out that the parametrization of the S-wave FSI there was not properly carried 
out because in the g 2 term in the amplitude, there are also D wave components. In 
Ref. [13], according to a unitarized chiral theory, an S wave effective Lagrangian and 
an effective scalar form factor were adopted. As a result, the nn invariant mass spectra 
of heavy quarkonium decays can be reproduced, but the cos 9* distribution still cannot 
be explained. 

In the T(3S) — > T(lS)n + n~ decay, the two-peak structure of the n + n~ invariant 
mass spectrum has been considered as a consequence of nn FSI [14, 15, 17, 18, 19] or 
the additional contribution from the D wave component of T(3S) [16]. Ignoring the 
contribution from the higher order pion momentum, Chakravarty et al. [17] explained 
the nn invariant mass spectrum of T(3S') — > T(lS')7r + 7r _ with x 2 /A r ( i./. = H-0/7 or 
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C.L.=14.0%, but only quantitatively gave the COS6 1 * distribution. Gallegos et al. [18] 
parametrized a more generalized amplitude in which both S and D wave contributions 
are included by fitting to the invariant mass distributions and the angular distributions 
of the decays mentioned at the beginning of this section. It is clear that the result of the 
parametrization should be consistent with the existing nn scattering data. Whether 
this condition has been satisfied in that calculation remains a question. Lahde et 
al. [19] showed that in order to better describe the ttti invariant mass spectra, the 
contribution of the pion re-scattering should be small in the ty(2S) — > J /^>-k + 'k~ or 
T(2S) -> T(lS)n + 7i- decay, but dominant in the T(3S) -> T(1S)ti + 7i- decay. Why 
it is so is still a puzzle. 

On the other hand, by introducing a bbqq resonance with J p = 1 + and mass of 
10.4-10.8 GeV, Anisovich et al. [20] explained the it + it~ invariant mass spectra of 
above mentioned bottomonium decays, but not the angular distribution of T(3S f ) — > 

T(15)7T + 7T-. 

To treat meson-meson S wave FSI properly, a chiral nonperturbative approach, 
called chiral unitary theory (ChUT), was recently proposed by Oiler and Oset [21] and 
later developed by themselves and others [22, 23, 24, 25, 26, 13, 27, 28, 29, 30, 31, 32] 
(for details, refer to the review article Ref. [33]). In this theory [21], the coupled channel 
Bethe-Salpeter equations (BSE), in which the lowest order amplitudes in the ChPT are 
employed as the kernel, are solved to resum the contributions from the s-channel loops 
of the re-scattering between pseudoscalar mesons. By properly choosing the three- 
momentum cutoff, the only free parameter in the theory, ChUT can well-describe the 
data of the S wave meson- meson interaction up to ^/s ~ 1.2GeV [21] which is much 
larger than the energy in the region where the standard ChPT is still valid, and can 
dynamically generate the scalar resonances a (/o(600)), /o(980) and ao(980) [21]. 

In this work, we adopt an amplitude used in Ref. [11], which includes both S and 
D wave contributions, and consider the S wave tttc FSI in the framework of ChUT 
to study the #(25) -> J/^tt+tt-, T(2S) -> T{lS)n + 7T-, T(3S) -> T(2S)ti + 7i- and 
T(3S) -> T(1S)ti + 7i- decays. 

The paper is organized as follows: In Section 2, the effective Lagrangian and ChUT 
are briefly introduced. In terms of the t-matrix written in Section 2, the ty(2S) — > 
J/^7y + 'k~ decay is discussed in Section 3. Section 4 is dedicated to the bottomonium 
7T + 7r~ transitions, and a brief summary is given in Section 5. 
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2 Brief formalism 

In the heavy quarkonium 7r + 7r~ transition process, the Lagrangian in the lowest 
order, which appropriately incorporates the chiral expansion with the heavy quark 
expansion, can be written as [11] 

£ = C + C S .B. (1) 

C = 9l A^B {v ^*Tr[{d v U){d p U) ] ] + g 2 A^B^*Tr[(v ■ dU)(v ■ dU) ] ] 

+ g z A^B^*Tr[{d^U){d u U) ] + {d^U) ] {d v U)\ + h.c. (2) 

Cs.b. = g±A^B {v) ^Tr[M{U + C/ f - 2)] 

+ ig,e^[v,A^d a B i ;> - (d,A^)v a B^>}Tr[M(U - £/+)] + h.c. (3) 

where gi denotes the coupling constants, U is a 3x3 matrix that contains the pseu- 
doscalar Goldstone fields, and M = diag{m u ,md,m s } is the quark mass matrix with 
m u , m d , m s being the masses of current quarks u, d and s, respectively. Affl and B^ 
are the fields of the initial and final states of heavy vector quarkonia, respectively, and 
v is the velocity vector of A. The tree diagram amplitude for the decay of a vector me- 
son into two pseudoscalar mesons and one vector meson in the rest frame of decaying 
particle can be expressed as [15, 17, 11] 

4 / / 

1 = —piidiPi ■ P2 + 92PiP° 2 + 93^1)6* ■ e + g^Vx^v + PivPi^e v \ (4) 

J 7T 

where f n = 93MeV is the decay constant of pion, p 1 and p 2 are the four-momenta of 
7r + and 7r~, respectively, p\ and p\ denote the energies of n + and ir~ in the lab frame, 
and e and e are the polarization vectors of the heavy quarkonia, respectively. It can be 
verified by the CLEO data [17] that by considering the chiral symmetry breaking scale 
and the heavy quark mass, the contribution from the last term (g 4 -term) is strongly 
suppressed [11]. Thus, the g 4 -term in Eq. (4) can be ignored and the amplitude can 
further be written as 

4 / 
V = —p(gm ■ P2 + 92p\pI + ff3"4K • £ ■ (5) 

It is noted that the D wave component exists in the gi term [10, 12]. Under Lorentz 
transformation, -p\ and p9, can be expressed as the functions of the momenta of pions 
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in the center of mass (cm.) frame of the nn system: 

rf = ^i^« + ^l|pt|cos^), (6) 

^ = ^=^(prH/3||Pi|cos0;), (7) 

where f3 is the velocity of the irir system in the rest frame of the initial particle, p\={pi* ■, 
p*) and P2 = (P2*> P2) are the four-momenta of n + and 7r~ in the cm. frame of the ixix 
system, respectively. So p\j>2 can t» e decomposed as 

P°iP°2 = I3£*?* 2 - ^Vo(cos^) - 2/3 2 pf P 2 (cos^)] (8) 

where P (cos9l) = 1 and P 2 (cos6 1 *) = |(cos 2 #* — |) are the Legendre functions of the 
0-th order and 2-nd order, respectively. 

Furthermore, the S wave nrr FSI which is important in this energy region should 
properly be included into the theoretical calculation. ChUT [21] is one of the suitable 
approaches for this job, because by using this theory, the S wave n — n scattering data 
up to 1.2GeV can be well reproduced. However, ChPT amplitudes in the 0(p 2 ) order 
are adopted as the kernel of the coupled-channel BSE [21], so the D wave FSI cannot 
be included. In the decays considered, the kinematical region is below 0.9GeV for the 
decay T(3S') — > T(1S)tt + 7t~ and below 0.6GeV for the others. One can see that they 
are far below the D wave resonant region (> 1.2GeV). So as a primary consideration, 
the D wave contribution comes only from the D wave terms appeared in Eq. (5). 

The basic diagrams for the V — > VPP decay, where V'(V) and P denote the vector 
and pseudoscalar mesons, respectively, are shown in Fig. 1. In the figure, (a) represents 
the V — > VPP decay without FSI, namely the tree diagram or Born term, and (b) 
describes the decay with tctt FSI. In Fig. 1(b), the -kit — > it + it~ t-matrix described 
by the solid black circle is obtained by the loop resummation [21], namely by a set 
of coupled-channel BSEs and both tttt and KK channels are included (for details, see 
Ref. [21]). To factorize the hit FSI from the direct V — > VPP decay part in Fig. 1(b), 
the on-shell approximation is adopted. Thus, only 7r + , 7r~, 7r° exist in the first loop 
which is directly linked to the V — > VPP vertex. The off-shell effects can be absorbed 
into the phenomeno logical coupling constants of the vertex. In fact, this approximation 
is often used in parameterizing the S wave FSI with phase shift data [17]. The full S 
wave 7T7T — > 7r + 7r~ t-matrix can be expressed as 

(tt+tt- + 7T-7T+ + TrVltlvr+TT-) = 2^=°_, (9) 
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FIG. 1: Diagrams of the vertex V'Vu + u~ without tttt FSI (a) and with tttt FSI (b). 



where t^ wn denotes the full S wave tttt — > tttt t-matrix in the isospin 1 = channel, 
which is the solution of a set of on-shell coupled channel BSEs [21]. Now, the t-matrix 
for the V' — > Vtt + tt~ decay can be expressed as 

t = V + V 0S -G-2t I =%, (10) 

where V is the amplitude of the tree diagram (a), V s is the S wave component of V 
and G is the two meson loop propagator 

G = i f d ' q 1 \ (ii) 

J (27r) 4 q 2 — ml + is (p — p — q) 2 — m 2 + is 

The numerical calculation is done by introducing a three-momentum cutoff q max - The 
value of the cutoff is taken as that used in [21], where the tttt scattering data can be 
well reproduced up to 1.2 GeV, namely, the tttt FSI in our model is consistent with 
the tttt scattering data. Moreover, this cutoff is also consistent with the dimensional 
regularization [22]. The analytic expression of the loop integral in Eq. (11) can be 
given as 

G = -i-{aarctan J- - ln[^(l + A)]}, (12) 
ott z Ac rrivr 



where a = J^-l and A = Jl + ^- 

* s y Qmax 



The differential decay width with respect to the tttt invariant mass and cos#* reads 

dm^dcosOl = 8M2(27r)3^ |t|2|pII|p31 (13) 

where J2 J2 describes the average over initial states and the sum over final states, and 
p 3 is the 3-momentum of the final vector meson in the lab frame. 
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3 Results for the tf(2S) -> J/^tt+tt" decay 

In the model, the parameters involved are the coupling constants g±, g 2 and g 3 . 
The values of the parameters can be determined by fitting the experimental data of 
the ty(2S) — > J/^!it + it~ process. It is shown that the resultant g% value is so small that 
we can safely take g 3 = 0. The remaining coupling constants g 1 and g 2 are obtained 
by fitting the total decay rate and the n + n~ invariant mass spectrum simultaneously 
The decay data of the V(2S) — > J/^n + n~ process are taken from ref.[34]. These BES 
data are normalized by using Y^^s) = 277keV and B(V(2S) -> J/^tt+tt") = 31.7% 
[35]. The resultant coupling constants are 

<?i = 0.106, g 2 /gi = -0.319, g 3 /gi = 0. (14) 



' iti>V '' / |t# ttH 

0.5 0.0 0.5 1.0 

cosej 

FIG. 2: ty(2S) — > J/^it + it~ decay, (a) it + it~ invariant mass spectrum, and (b) cos#* 
distribution with 9* being the angle between the moving direction of n + in the 7i + n~ 
cm. frame and the moving direction of the 7r + 7r~ system in the lab system. The data 
points are taken from the BES result [34]. The solid curves denote our best fitted 
results. 

Our best fit to the 7r + 7r~ invariant mass spectrum and the cos#* distribution are 
shown in Fig. 2. It is found that the ii + 7i~ invariant mass spectrum is well-fitted, but 
the theoretical angular distribution is somewhat too flat. Note that in fitting angular 
distribution, we only consider cos#* from -0.8 to 0.8, because the efficiency correction 
to the data at large |cos6>*| is not accurate enough [34]. The deviation in angular 
distribution implies that the D wave contribution is somehow too small. In fact, as 
discussed in Section 2, in our calculation, the D wave FSI is not included. It is also 
found that the S wave FSI enhances the invariant mass spectrum considerably. It 
should be noted that due to P 2 (cos 9*) = 1/2 (cos 2 0* — 1/3), integrating over cos 0* will 
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m (GeV) 



result that the D wave contribution is not so important in the invariant mass spectrum. 
However, in the angular distribution, the cos ^-dependence, and consequently the D 
wave effect, will explicitly show up. Thus, we deem that the deviation in angular 
distribution may be due to lack of D wave FSI. In fact, it can be confirmed in the 
following way: In Ref.[34, 11], without FSI, the authors can well- reproduce both the 
7T7T invariant mass spectrum and the angular distribution with gi = 0.30±0.01, g2/gi= 
—0.35 ± 0.03 and g% = 0. However, with the S wave FSI, in the best data fitting to 
the 7T7T invariant mass spectrum, the resultant g\ is 0.106 which is 3 times less than 
that in Ref. [11], while g^j Q\ keeps almost the same value as that in Ref. [11]. This 
means that the effect of the S wave FSI is so large that g\ has to be much smaller to 
explain the mi invariant mass spectrum. As a consequence, the D wave component is 
also greatly reduced. If we naively multiply a factor of 3 to the D wave component, 
the angular distribution can also be reproduced better. Unfortunately, as mentioned 
above, the D wave FSI can not be treated properly in the simple ChUT approach [21]. 

4 The bottomonia tt + tt~ transitions 

Similar to the case of V(2S) -> J/^ti + ti~, in the T(2S) -> T (I^tt+tt- or T(3S) -> 
T(2S)it + it~ decay, g 3 = can be adopted. But for the T(3S) — > Y(IS)it + it~ decay, 
the S wave FSI is no longer a main contributor, and a finite value of gs is requested. 
With this consideration, Mannel et al. showed that the resultant values of gil Q\ for the 
T(2S) — > Y(IS)it + it~ and T(3S) — > Y(2S)it + it~ processes are very close, but quite 
different from that for the T(3S) — > T(1S)tt + 7t~ process. The latter one is about ten 
times larger than the former [11]. This is somewhat unnatural. Suffice to say, the pions 
involved in T(3S) — > T(lS)n + ir~ are somewhat harder than in the other bottomonium 
transitions, and in principle the values of ^/^i f° r these processes should not be the 
same for different dynamical regions. However, in the decay processes considered, the 
vector mesons involved are all in the S wave state, and the particles involved are in the 
same mass scale, and the difference among kinematical regions is not too large. Thus 
we deem that the values of g^l Q\ for these processes should be very close. To reduce the 
number of free parameters, we take the same g<ijg\ value for different T(nS) decays. 

The decay data for T(2S) — > T(1S)tt + 7t~ are taken from [36] and for T(3S) — > 
T(2S)tc + 7t~ and T(3S) — > T(lS')7r + 7r _ from [1]. To get the physical coupling con- 
stants, the data for T(2S') — > T(lS')7r + 7r _ [36] is normalized by r X (2S) = 43/ceV and 
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FIG. 3: (a) The ir + ir~ invariant mass spectrum for the T(2S) — > T(lS)ir + ir~ process, 
(b) the cos 9* distribution for the T(2S) — > Y(lS')7r + 7r _ process, and (c) the 7r + 7r~ 
invariant mass spectrum for the T(3S) — > T(2S')7r + 7r~ process. 



B(T(2S) -> T(l 1 S)7T+7r-) = 19.2% [35]. Our calculated results are plotted in Figs. 3-4. 
It is shown that the resultant ittt invariant mass spectra for both T(2S) — > T(lS')7r + 7r _ 
and T(3S) — > T(2S')7r + 7r _ decays agree with the data values, but the angular distribu- 
tion for the former one is somewhat flat, which might also be due to the same reason 
discussed in Section 3. On the other hand, there is almost no way to fit the 7r + 7r~ 
invariant mass spectrum and the cos#* distribution of the T(3S) — > T(lS)ir + ir~ pro- 
cess simultaneously, even if gijg\ is further released as a free parameter. The resultant 
parameters are listed in Tab. 1. 



TABLE 1: Resultant parameters through data fitting. 



Decay 


9i 


92/91 


93/91 


T(2S) -> T{1S)ti + 7i- 


0.0944 


-0.230 





T(3S) -> T{2S)tt + 7i- 


0.768 


-0.230 





T{3S) -> T(1 1 S)7T + 7T- 


0.0123 


0.564 


-13.602 
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FIG. 4: (a) The n + ir~ invariant mass spectrum and (b) the cos6>* distribution for the 
T(3S) -> T(1S)ti + 7i- process. 

4.1 Sequential decay mechanism 

In order to explain the decay data of the T(3S) — > T(IS)it + it~ decay, we propose 
an additional sequential decay mechanism where an intermediate state, called X, is 
introduced. Additional Feynman diagrams for T(nS) — > T(mS)7r + 7i^ are shown in 
Fig. 5, where (a) depicts the tree diagram and (b) the diagram including the nn S wave 
FSI. We adopt a simple S wave coupling for T(nS) — > nX. The quantum numbers of 
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FIG. 5: Diagrams of the sequential X decay mechanism (a) at tree level, and (b) with 
the tttt FSI. 



X should be J p = 1 + and 1 = 1. The decay amplitude of Fig. 5(a) can be written as 

v tree „ • * ( + f X+ Px + /™X + P X -P X - / m X , 

" p x+ ~ m x + imxLx P x - ~ m x + i™> x r x 
where px+ and p x - are the momenta of X + and X~ respectively, g nm is an effective 
coupling constant among T(nS), T(mS), n + and 7T~ via an intermediate resonant 
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state X. In fact, g nm is the product of two coupling constants g n x and g m x where 
9kx(k = n,m) denotes the coupling constant for the T(kS)Xir vertex. To further 
consider the effect of the irn S wave FSI, the contribution of Fig. 5(b) should be 
included. In this figure, the three-propagator loop can be expressed as 

q\xv , f dq* -g^+ PxPx /m x 1 _ 1 

x J (27r) 4 p 2 x — m x + ie q 2 — ml + ie (p — p — q) 2 — ml + ie 

where px = P — q is the four- momentum of X. The calculation is carried out in the 
cm. frame of the 7T7T system with the same cutoff value used in the two-meson loop 
calculation. As argued in ref.[20], terms with e'^tPxPx/ m x m Eqs. 15 and 16 can 
be neglected, because of the expected heavy mass of X. Then the total t-matrix can 
finally be written as 

t = V + V 0S ■ G ■ 2t[Z% + + 9nJ/ v G^ ■ K=%. (17) 

4.2 Results for bottomonium 7r + 7r transitions 

In terms of the t-matrix in Eq. (17), we calculate the 7r + 7r~ invariant mass spectra 
and the cos^* distributions of the bottomonium 7r + 7r~ transitions. Similar to the 
argument given in the former sections, we take g 3 = for the T(2S) — > Y(IS)it + it~ and 
T(3S) — > T(2S')7r + 7r _ decays and g 3 as a free parameter for the T(3S) — > T(lS')7r + 7r _ 
process. We also demand the values of (^/^i to be the same for all three decays for 
reducing the number of free parameters. The values of g\ and g 2 are determined by 
fitting the experimental decay data [36, 1]. 

The calculated results are plotted in Fig. 6. It is shown that not only both the 
7T7T invariant mass spectrum and the cos#* distribution of the T(3S) — > T(lS')7r + 7r _ 
process can simultaneously be well explained, but also a consistent description of other 
bottomonium 7r + 7r~ transitions can be obtained. The resultant parameters are tabu- 
lated in Tab. 2. 

To understand thoroughly the roles of different terms in the 7r + 7r~ invariant mass 
spectrum and the cos#* distribution of the T(3S) — > Y(lS)ir + ir~ decay, it is neces- 
sary to analyze their individual contributions. The results are shown in Fig. 7. In the 
figure, the solid curves represent our best fitted results, and the dotted, dashed, and 
dash-dotted curves describe the contributions from the terms without X and with X 
only and the interference term respectively, and the dash-dot-dotted curves represent 
the tree level contributions with X only. The calculated ttti invariant mass spectrum 
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FIG. 6: The 1111 invariant mass spectra and the cos6>* distributions in the T(nS) 
Y(mS)n + iT~ decays. 



(Fig. 7 (a)) shows that the contribution from X plays a dominant role, the contri- 
bution from the terms without X can qualitatively but not quantitatively give the 
two-peak feature, and the interference term contributes constructively in the smaller 
tttc invariant mass region but destructively in the larger invariant mass region. The 
resultant cos#* distribution (Fig. 7 (b)) further shows that the contribution from X, 
even in the tree level, produces almost the whole angular distribution structure. Al- 
though the scalar meson a dynamically generated by the S wave tttc FSI in ChUT 
[21] can make a peak around its pole position at about 450 MeV in the 7T7T invariant 
mass spectrum, the contribution from the diagrams without X is not dominant due to 
smaller values of coupling constants. Thus, an additional D wave FSI which provides 
a flat contribution in the invariant mass region considered will not be an important 
contributor. These indicate that the intermediate state X is very important in repro- 
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TABLE 2: Resultant parameters in the data fitting. In the first column, n — > m 
denotes the T(nS) — > Y(mS)ir + 7r~ decay. 



Decay 


0i 


92/91 


93/91 


<w(GeV 2 ) 


m x (GeV) 


Tx(GeV) 


2 -> 1 


0.0886 


-0.230 





-2.316 


10.080 


0.655 


3^2 


0.769 


-0.230 





-0.00418 


3 -> 1 


0.00546 


-0.230 


4.949 


4.712 




0.4 0.6 ' 0.8 ~ " ' 0.0 " ~ 0.5 * " 1.0 



(GeV) cose " 

FIG. 7: The contributions from different components to (a) the n + n~ invariant mass 
spectrum and (b) the cos#* distribution in the T(3S) — > Y(lS)ir + ir~ decay. The solid 
curves represent our best fitted results, the dotted, dashed and dash-dotted curves 
describe the contributions from terms without X, with X only and the interference 
term, respectively. The dash-dot-dotted curves represent the tree level contributions 
with X only. 

ducing not only the 7r + 7r~ invariant mass spectrum but also the cos #* distribution in 
the T(3S) -> T(1S)ti + 7i- decay. 

If we further consider the quark structures of the particles involved, the intermediate 
state should contain b, b, q and q. This state might be a tetraquark state, for instance, 
bbud, with J p = 1 + and I — 1 for X + or a BB bound state, for instance, B + B°, for 
X+. 

It should be mentioned that similar mechanism was also proposed by V. V. Anisovich 
et al. [20]. In their paper, a trivial S wave coupling was used in the effective vertex 
T (nS)T \mS)ixn which is described by Eq. (5) in our model with both S wave and D 
wave components. With that mechanism, they successfully reproduced the nn invariant 
mass spectra of the T(nS) — > T(mS')7r + 7r _ decays, but did not give reasonable cos#* 
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FIG. 8: Dalitz plot for T(3s) — > T(ls)7r + 7r . The mass squared of T(1S)tt system is 
plotted against the mass squared of the n + n~ system. 



distributions due to the dominance of the tttt S wave in their model. As a result, the 
estimated mass of the additional intermediate state is in the range of 10.4-10.8 GeV 
which is located outside the data area of the Dalitz plot where the data points in the 
direction of m X7r are located from 9.6 GeV to 10.2 GeV. Thus the effect of the state 
does not show up in the Dalitz plot. 

The mass and width of the intermediate state in our work are different from those 
in Ref. [20]. We also present the Dalitz plot for the T(3S) — > T(lS')7r + 7r _ decay in 
Fig. 8. It is shown that although the estimated mass of X in our model (M x = 10.08 
GeV) is inside the data area in the Dalitz plot, the signal of X in the Dalitz plot is 
not very clear due to its large width of 0.655 GeV. It does not conflict with the CLEO 
experiment [36]. Moreover, we would mention that with the typical values, Mx=10.5 
GeV and Fx=0.15 GeV given in Ref. [20], we cannot produce a cos#* distribution that 
is consistent with the experimental data [36]. 

5 Summary 

Starting from an effective Lagrangian and further employing ChUT to include 
the 7T7T S wave FSI properly, the tt + tt" transitions of heavy quarkonia are inten- 
sively studied. In order to consistently explain the 7r + 7r~ invariant mass spectra 
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and angular distributions in the mentioned processes simultaneously, especially in the 
T(3S) — > Y(lS)n + 7i~ decay process, an additional sequential process, where an inter- 
mediate state X is introduced, is further considered in the bottomonium transitions. 
With such a process included, all the 7r + 7r~ transition data can be well-explained, es- 
pecially the angular distribution of the T(3S) — > T(lS)ir + ir~ decay. As a consequence, 
the newly introduced intermediate state should have quantum numbers of J p = 1 + and 
I — 1, a mass of about 10.08 GeV and a width of about 0.655 GeV. The quark content 
of the state should be bbqq. It might be a tetraquark state or a BB bound state. The 
detailed inner structure of the state should be carefully studied both theoretically and 
experimentally. 
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